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j3 ' This paper gives an overview of recent results concerning the long time dy- 

namics of repeated interaction quantum systems in a deterministic and ran- 
dom framework. We describe the non equilibrium steady states (NESS) such 
systems display and we present, as a macroscopic consequence, a second law 
of thermodynamics these NESS give rise to. We also explain in some details 
the analysis of products of certain random matrices underlying this dynamical 
problem. 
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1. Introduction and Model 

00 

A repeated interaction quantum system consists of a reference quantum 
subsystem S which interacts successively with the elements £ m of a chain 
C = £ i +£2 + - ■ ■ of independent quantum systems. At each moment in time, 
. S interacts precisely with one £ m (m increases as time does) , while the other 

elements in the chain evolve freely according to their intrinsic dynamics. 
The complete evolution is described by the intrinsic dynamics of S and of 
all the £ m , plus an interaction between S and £ m , for each m. The latter is 
characterized by an interaction time r m > 0, and an interaction operator V m 
(acting on S and £ m )\ during the time interval [t\ + ■ ■ •+r m _i, t\ + ■ ■ ■+r m ), 
iS is coupled to £ m only via V rn . Systems with this structure are important 
from a physical point of view, since they arise naturally as models for funda- 
mental experiments on the interaction of matter with quantized radiation. 
As an example, the "One atom maser" provides an experimental setup in 
which the system S represents a mode of the electromagnetic field, whereas 
the elements £k describe atoms injected in the cavity one by one, which 
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interact with the field during their flight in the cavity. After they leave the 
cavity, the atoms encode some properties of the field which can be mea- 
sured on these atoms 14 , 16 For repeated interaction systems considered as 
ideal, i.e. such that all atoms are identical with identical interactions and 
times of flight through the cavity corresponding mathematical analyses are 
provided in 17 , To take into account the unavoidable fluctuations in the 
experiment setup used to study these repeated interaction systems, mod- 
elizations incorporating randomness have been proposed and studied in 8 
and. 9 With a different perspective, repeated quantum interaction models 
also appear naturally in the mathematical study of modelization of open 
quantum systems by means of quantum noises, see 4 and references therein. 
Any (continuous) master equation governing the dynamics of states on a 
system S can be viewed as the projection of a unitary evolution driving the 
system S and a field of quantum noises in interaction. It is shown in 4 how 
to recover such continuous models as some delicate limit of a discretization 
given by a repeated quantum interaction model. Let us finally mention 15 
for results of a similar flavour in a somewhat different framework. 

Our goal is to present the results of the papers 7 , 8 and 9 on (random) 
repeated interaction quantum systems, which focus on the long time be- 
haviour of these systems. 

Let us describe the mathematical framework used to describe these 
quantum dynamical systems. According to the fundamental principles of 
quantum mechanics, states of the systems S and £ m are given by normal- 
ized vectors (or density matrices) on Hilbcrt spaces TLs and respec- 
tively, 3 ' 621 . We assume that dimWs < oo, while dimHg^ may be infinite. 
Observables Ag and Ag m of the systems S and £ m are bounded opera- 
tors forming von Neumann algebras Wis C B(Hs) and Wls m C B(Hs m ). 
They evolve according to the Heisenberg dynamics Isin ag(As) and 
M.3 t a £ (Ag m ) , where a l s and a £ are *-automorphism groups of 9Jts 
and 9Jl£ m , respectively, see e.g. 6 We now introduce distinguished reference 
states, given by vectors ips 6 and ^ps m £ ^£ m - Typical choices for ips, 
ip£ m are equilibrium (KMS) states for the dynamics a s , ct £m , at inverse 
temperatures /?s, 0£ m . The Hilbert space of states of the total system is 
the tensor product 

H = H S ®H C , 



a A normalized vector tp defines a "pure" state A i— ► (tp,Aip) = Tr(g^A), where q^, = 
|i/>)(i/i|. A general "mixed" state is given by a density matrix q = Y2 n >lP n @i>n> where 
the probabilities p n > sum up to one, and where the tp n are normalized vectors. 
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where Tic = ® m >i^ m i an d where the infinite product is taken with 
respect to Vc = ® m >i^p£ m - The non-interacting dynamics is the prod- 
uct of the individual dynamics, defined on the algebra Wig ® m>1 %ft£ m 
by a s®m>i a £ m - I* wm P rove useful to consider the dynamics in the 
Schrodinger picture, i.e. as acting on vectors in 7i. To do this, we first 
implement the dynamics via unitaries, satisfying 

a* # (A#) = e itL *A # e~ HL *, t G K, andL # V># = 0, (1) 

for any A# G where # stands for either S or £ m . The self-adjoint 

operators Ls and Ls m , called Liouville operators or "positive temperature 
Hamiltonians" , act on TLs and He m , respectively. The existence and unique- 
ness of L# satisfying (1) is well known, under general assumptions on the 
reference states V'#- 6 We require these states to be cyclic and separating. 
In particular, (1) holds if the reference states are equilibrium states. Let 
r m > and V m G 9Jts <g> %Jls m be the interaction time and interaction oper- 
ator associated to S and £ m . We define the (discrete) repeated interaction 
Schrodinger dynamics of a state vector <p G H, for m > 0, by 

U(m)<f> = c- iZ ™ ■ ■ ■ c- ii2 c- [il (f>, (2) 

where 

Lk = i~kL k + T fe Ls n (3) 

describes the dynamics of the system during the time interval [t\ + • • • + 
Tk-i,T\ + • • • + Tfe), which corresponds to the time-step k of our discrete 
process. Hence Lk is 

L k = L s + L £k +Vk, (4) 

acting on Hs <8> ~He k ■ We understand that the operator Lg n in (3) acts 
nontrivially only on the n-th factor of the Hilbcrt space Tic of the chain. 
As a general rule, we will ignore tensor products with the identity operator 
in the notation. 

A state £?(•) = Tr(p • ) given by density matrix p on H is called a normal 
state. Our goal is to understand the large-time asymptotics (to — ► oo) of 
expectations 

q (U(m)*OU(m)) = g{a m {0)), (5) 

for normal states g and certain classes of observables O that we specify 
below. We denote the (random) repeated interaction dynamics by 

a m (0) = U{m)*OU{m). (6) 
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1.1. Van Hove Limit Type Results 

A first step towards understanding the dynamics of repeated interaction 
quantum systems reduced to the reference system S was performed in 
the work. 2 This paper considers Ideal Repeated Quantum Interaction Sys- 
tems which are characterized by identical elements = £ in the chain 
C, constant interaction times Tk = t and identical interaction operators 
Vfe = V E Wis ® 9tt<£ between S and the elements £ of the chain. In 
this setup, a Van Hove type analysis of the system is presented, in sev- 
eral regimes, to describe the dynamics of observables on S in terms of a 
Markovian evolution equation of Lindblad type. Informally, the simplest 
result of 2 reads as follows. Assume the interaction operator V is replaced 
by XV, where A > is a coupling constant, and let m, the number of inter- 
actions during the time T = mr, scale like m c± t/X 2 , where < t < oo and 
r are fixed. Assume all elements of the chain are in a same thermal state 
at temperature (3. Then, the weak coupling limit A — > of the evolution of 
any observable O acting on S obtained by tracing out the chain degrees of 
freedom from the evolution (6) satisfies, after removing a trivial free evolu- 
tion, a continuous Lindblad type evolution equation in t. The temperature 
dependent generator is explicitely obtained from the interaction operator 
V and the free dynamics. The asymptotic regimes of the parameters (A, t) 
characterized by r — > and tA 2 < 1 are also covered in, 2 giving rise to dif- 
ferent Lindblad generators which all commute with the free Hamiltonian on 
S. The critical situation, where r — > with rA 2 = 1 yields a quite general 
Lindblad generator, without any specific symmetry. In particular, it shows 
that any master equation driven by Lindblad operator, under reasonable 
assumptions, can be viewed as a Van Hove type limit of a certain explicit 
repeated interaction quantum system. 

By contrast, the long time limit results obtained in 7 , 8 and 9 that we 
present here are obtained without rescaling any coupling constant or pa- 
rameter, as is usually the case with master equation techniques. It is pos- 
sible to do without these approximations, making use of the structure of 
repeated interaction systems only, as we now show. 

2. Reduction to Products of Matrices 

We first link the study of the dynamics to that of a product of reduced dy- 
namics operators. In order to make the argument clearer, we only consider 
the expectation of an observable A$ e Wig , and we take the initial state of 
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the entire system to be given by the vector 

ipo = ips ®V>c, (7) 

where the ips and ipc are the reference states introduced above. We'll com- 
ment on the general case below. The expectation of A$ at the time-step m 
is 

(Vo, °r(A s )4, Q ) = (Vo, Pc iZl ■ ■ ■ e iZ ™A s c- iZ ™ • • • e- i£l P^ ) , (8) 
where we introduced 

P = l ns <g)P^ m , (9) 

m>l 

the orthogonal projection onto 7Ys<S>CV>c- A first important ingredient in our 
analysis is the use of C-Liouvilleans introduced in 11 , which are operators 
Kk defined by the properties 

e iZk Ac- iik = e iKk Ae~ iKk , (10) 

K k 4> s ® i>c = 0, (11) 

where A in (10) is any observable of the total system. The identity (10) 
means that the operators Kk implement the same dynamics as the Lk 
whereas relation (11) selects a unique generator of the dynamics among 
all operators which satisfy (10). The existence of operators Kk satisfying 
(10) and (11) is rooted to the Tomita-Takesaki theory of von Neumann 
algebras, c.f. 11 and references therein. It turns out that the Kk are non- 
normal operators on H, while the Lk are self-adjoint. Combining (10) with 
(8) we can write 

(^a m (A s )^) = (iPo,Pe iK ^--e iK ™PA s ip ). (12) 

A second important ingredient of our approach is to realize that the inde- 
pendence of the sub-systems £ m implies the relation 

Identifying Pc lKk P with an operator on Hs, wc thus obtain from (12) 
and (13), 

(Vo, a m (A s )^ ) — (ips, Mi ■ ■ ■ M m A s fo) . (14) 

It follows from (11) that Mktps = 4>s, f° r all fc, and, because the operators 
Mk = Pe lKk P implement a unitary dynamics, we show (Lemma 4.1) that 
the Mk are always contractions for some suitable norm || • ||| on C d . This 
motivates the following 



June 1, 2008 15:8 



WSPC - Proceedings Trim Size: 9in x 6in qmathlO 



6 

Definition: Given a vector ijjg G C d and a norm on ||| • ||| on C d , we call 
Reduced Dynamics Operator any matrix which is a contraction for ||| • || 
and leaves ips invariant. 

Remark: In case all couplings between S and £ k are absent, V k = 0, 
Mfe = e lTkLs is unitary and admits 1 as a degenerate eigenvalue. 

We will come back on the properties of reduced dynamics operators 
(RDO's, for short) below. Let us emphasize here that the reduction process 
to product of RDO's is free from any approximation, so that the set of 
matrices {Mk = Pe lKk P}keN encodes the complete dynamics. In particu- 
lar, we show, using the cyclicity and separability of the reference vectors 
tps, tp£ k , that the evolution of any normal state, not only (ipQ, ■ ipo), can be 
understood completely in terms of the product of these RDO's. 

We are now in a position to state our main results concerning the asymp- 
totic dynamics of normal states g acting on certain observables O. These 
result involve a spectral hypothesis which we introduce in the next 

Definition: Let M(e) denote the set of reduced dynamics operators whose 
spectrum a{M) satisfies a(M) D {z & C \ \z\ = 1} = {1} and 1 is simple 
eigenvalue. 

We shall denote by P\.m the spectral projector of a matrix M corre- 
sponding to the eigenvalue 1. As usual, if the eigenvalue 1 is simple, with 
corresponding normalized eigenvector ips, we shall write P\,m = |^s)(V'l 
for some ip s.t. (ip\ips) = 1- 

3. Results 

3.1. Ideal Repeated Interaction Quantum System. 

We consider first the case of Ideal Repeated Interaction Quantum Systems, 
characterized by 

£ k = £, L £k = I £ , V k = V, T k =r for all k > 1, 
M k = M, Vk> 1. (15) 

Theorem 3.1. let a n be the repeated interaction dynamics determined by 
one RDO M. Suppose that M G M( E ) so that P^.m = \i>s)(ip\- Then, for 
any < 7 < inf zecr( - M )\{ 1 }(l — any normal state g, and any As G Wis, 

g (a n (A s )) - (V, Asi>s) + 0(e"^). (16) 
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Remarks: 1. The asymptotic state {tp\ ■ tps) and the exponential decay rate 
7 are both determined by the spectral properties of the RDO M. 

2. On concrete examples, the verification of the spectral assumption on 
M can be done by rigorous perturbation theory, see. 7 It is reminiscent of a 
Fermi Golden Rule type condition on the efficiency of the coupling V, see 
the remark following the dcfiniton of RDO's. 

3. Other properties of ideal repeated interaction quantum systems are 
discussed in 7 , e.g. continuous time evolution and correlations. 

For deterministic systems which are not ideal, the quantity g(a n (As)) 
keeps fluctuating as n increases, which, in general, forbids convergence, see 
Proposition 5.3. That's why we resort to ergodic limits in a random setup, 
as we now explain. 

3.2. Random Repeated Interaction Quantum System 

To allow a description of the effects of fluctuations on the dynamics of 
repeated interaction quantum systems, we consider the following setup. 

Let uj i— > M(uj) be a random matrix valued variable on C d defined on a 
probability space (O, J 7 , p). We say that M(u) is a random reduced dynamics 
operator (RRDO) if 

(i) There exists a norm || • ||| on C d such that, for all ui, M(ui) is a 
contraction on C d for the norm 1 1 1 • 1 1 . 

(ii) There exists a vector ^5, constant in ui, such that M(u))ips — *Ps, 
for all uj. 

To an RRDO ui 1— ► M(cu) on Q is naturally associated a iid random 
reduced dynamics process (RRDP) 

aJi-»M(wi)---M(w n ), Zuen w \ (17) 

where we define in a standard fashion a probability measure dP on f2 N by 

dP = n^idpj, where d Pj - = dp, Vj E N*. 

We shall write the expectation of any random variable / as E[/]. 

Let us denote by oQ, ZJ € £1 N , the process obtained from (6), (14), 
where the Mj = M(ujj) in (14) are iid random matrices. We call the 
random repeated interaction dynamics determined by the RRDO M(w) = 
Pe lK ^P. It is the independence of the successive elements £k of the chain 
C which motivates the assumption that the process (17) be iid. 
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Theorem 3.2. Let be the random repeated interaction dynamics de- 
termined by an RRDO M(lu). Suppose that p(M(w) G M(e)) > 0. Then 
there exists a set ft C N , s.t. P(O) = 1, and s.t. for any uJ e 0, any 
normal state g and any As E Wis, 

1 N 

^jv£<?(<^)) = W^s>, (18) 
where 9 = P* >m[M] i>s- 

Remarks: 1. Our setup allows us to treat systems having various sources 
of randomness. For example, random interactions or times of interactions, 
as well as random characteristics of the systems £ m and S such as random 
temperatures and dimensions of the £ m and of S. 

2. The asymptotic state (9, ■ ips) is again determined by the spectral 
data of a matrix, the expectation E[M] of the RRDO M(u). Actually, our 
hypotheses imply that E[M] belongs to M(e), see below. 

3. The explicit computation of the asymptotic state, in this Theorem 
and in the previous one, is in general difficult. Nevertheless, they can be 
reached by rigorous perturbation theory, see the examples in 7 , 8 and. 9 

3.3. Instantaneous Observables 

There are important physical observables that describe exchange processes 
between S and the chain C and, which, therefore, are not represented by 
operators that act just on S. To take into account such phenomena, we 
consider the set of observables defined as follows. 

Definition: The instantaneous observables of S + C are of the form 

= A s 2> r j= _ l BV>, (19) 
where A s £ Ms and G M £m+] . 

Instantaneous observables can be viewed as a train of / + r + 1 observables, 
roughly centered at £ m , which travel along the chain C with time. 

Following the same steps as in Section 2, we arrive at the following 
expression for the evolution of the state ipo acting on an instantaneous 
observable O at time m: 

(Vo, a m {0)i> Q ) = (Vo, PM l ■ ■ ■ M m _,_iJV m (0)P^o) . (20) 
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Here again, P is the orthogonal projection onto Tts, along ipc- The operator 
N m {0) acts on Hs and has the expression (Proposition 2.4 in 9 ) 

N m (0)^ Q = (21) 

p e ir m _,Z m _, . . .jr m L m ( As r = _ ; B U)^-ir m L m . . . e -ir m _,Z m _ ^ 

We want to analyze the asymptotics m — > oo of (20), allowing for ran- 
domness in the system. We make the following assumptions on the random 
instantaneous observable: 

(Rl) The operators Mk are RRDO's, and we write the corresponding iid 
random matrices Mk = M(uik), fc = 1, 2, • ■ • , . 

(R2) The random operator N m (0) is independent of the Mk with 1 < 
k < m — I — 1, and it has the form N(oj m -i 7 . . . , u m+r ), where 
N : f2 r+i+1 — > B(C d ) is an operator valued random variable. 

The operator Mk describes the effect of the random fc-th interaction on 
5, as before. The random variable N in (R2) does not depend on the time 
step m, which is a condition on the observables. It means that the nature 
of the quantities measured at time m are the same. For instance, the Bm 
in (19) can represent the energy of £ m +j, or the part of the interaction 
energy V m +i belonging to £ m +j, etc. Both assumptions are verified in a 
wide variety of physical systems: we may take random interaction times 
Tfc = r(u)k), random coupling operators T4 = V(uik), random energy levels 
of the £k encoded in Lg k = Ls(u)k), random temperatures f3g k = Ps(^k) of 
the initial states of £ k , and so on. 

The expectation value in any normal state of such instantaneous ob- 
servables reaches an asymptotic value in the ergodic limit given in the next 

Theorem 3.3. Suppose that p(M(u) E M-(e)) 0- There exists a set 
Q C il N of probability one s.t. for any w e O, for any instantaneous 
observable O, (19), and for any normal initial state g, we have 

lim - V e«(0)) - (0,E[iV]V>s> , nN} e Wl s . (22) 

fi^oo n *■ — ' 
m=l 

Remarks 1. The asymptotic state in which one computes the expectation 
(w.r.t the randomness) of N is the same as in Theorem 3.2, with 8 = 

-f > l*E[M]V'5- 

2. In case the system is deterministic and ideal, the same result holds, 
dropping the expectation on the randomness and taking 9 = tp, as in The- 
orem 3.1, see. 7 
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3.4. Energy and Entropy Fluxes 

Let us consider some macroscopic properties of the asymptotic state. The 
systems we consider may contain randomness, but we drop the variable ZU 
from the notation. 

Since we deal with open systems, we cannot speak about its total en- 
ergy; however, variations in total energy are often well defined. Using an 
argument of one gets a formal expression for the total energy which is 
constant during all time-intervals [t to _i,t to ), and which undergoes a jump 

j(m) :=a m (V m+1 -V m ) (23) 

at time step m. Hence, the variation of the total energy between the instants 
and m is then AE(m) = J^fcLi iW- The relative entropy of g with respect 
to g , two normal states on Wl, is denoted by Ent(g\g ). Our definition 
of relative entropy differs from that given in 6 by a sign, so that in our 
case, Ent(g\g ) > 0. For a thermodynamic interpretation of entropy and 
its relation to energy, we assume for the next result that ips is a (pts,a s )- 
KMS state on 9Jt 5 , and that the tp £m are (f3 £fn , a| m )-KMS state on M £m , 
where [3s is the inverse temperature of S, and [3 £m are random inverse 
temperatures of the £ m . Let go be the state on Wl determined by the vector 
V'o = ips ® tpc — ips ® m ip£ m ■ The change of relative entropy is denoted 
AS(m) := Ent(f) o a m \g ) — Ent(p|p )- This quantity can be expressed in 
terms of the Liouvillean and interaction operators by means of a formula 
proven in. 12 

One checks that both the energy variation and the entropy variations 
can be expressed as instantaneous observables, to which we can apply the 
results of the previous Section. Defining the asymptotic energy and entropy 
productions by the limits, if they exist, 

/AE{m\ 1Tn , AS(m) 
lim q =: AE+ and lim =: dS+, (24) 

m— >co \ m J m^oo m 

we obtain 

Theorem 3.4 (2 nd law of thermodynamics). Let g be a normal state 
on Wl. Then 

dE+ = (6,E[P(L S + V - e iTL (L s + V)e-' 1TL )P]ip s ) a.s. 
dS+ = (e,E[0 e P(Ls +V- e iTL (L s + V)c- WL )P]^ S ) a.s. 

The energy- and entropy productions dE + and d5+ are independent of the 
initial state g. If f3 £ is deterministic, i.e., uj -independent, then the system 
satisfies the second law of thermodynamics: dS+ — (3 £ dE + . 
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Remark: There are explicit examples in which the entropy production can 
be obtained via rigorous perturbation theory and is proven to be strictly 
positive, a sure sign that the asymptotic state is a NESS, see 7 . 

As motivated by (14), the theorems presented in this Sections all rely 
on the analysis of products of large numbers of (random) RDO's. The rest 
of this note is devoted to a presentation of some of the key features such 
products have. 

4. Basic Properties of RDO's 

Let us start with a result proven in 7 as Proposition 2.1. 

Lemma 4.1. Under our general assumptions, the set of matrices {Mj}j e ^» 
defined in (14) satisfy Mjips — ips, for all j E N*. Moreover, to any <f> E TLs 
there corresponds a unique A E 9Jl$ such that <f> = Alps- \\\4>\\\ '■— ll-^-lle(w s ) 
defines a norm on Hs, and as operators on Hs endowed with this norm, 
the Mj are contractions for any j E W . 

Again, the fact that ips is invariant under Mj is a consequence of (11) and 
their being contractions comes from the unitarity of the quantum evolution 
together with the finite dimension of Tis- 

As a consequence of the equivalence of the norms || • || and ||| • |||, we get 

Corollary 4.1. We have 1 E a(Mj) C {z | \z\ = 1} and 

sup {|| M jn M 3n _ 1 ■■■M jl \\, n E N*, j k E N*} = C < oo 

Actually, if a set of operators satisfies the bound of the Corollary, it 
is always possible to construct a norm on C d relative to which they are 
contractions, as proven in the next 

Lemma 4.2. Let R = {Mj E Md(C)}j e j, where J is any set of indices 
and C(R) > 1 such that 

\\M n M n ■ ■ ■ M jn || < C(R), V{ji}i=i,... , n G J n , Vn E N. (25) 

Then, there exists a norm ||| • ||| on C d , which depends on R, relative to 
which the elements of R are contractions. 

Proof: Let us define T C M d (C) by 



T = U„ eN U 0l;j2; ... jn)eJ „ M n M n ■ ■ ■ Mj, 



(26) 



June 1, 2008 15:8 



WSPC - Proceedings Trim Size: 9in x 6in qmathlO 



12 



Obviously R C T, but the identity matrix I does not necessarily belong to 
T. Moreover, the estimate (25) still holds if the Mj i 's belong to T instead 
of R. For any ip e C d we set 

HMH = sup \\M<p\\ > (27) 
m eTui 

which defines a bona fide norm. Then, for any vector ip and any element N 
of T we compute 

HliV^III = sup \\MN<p\\< sup \\M<p\\ = |IMII> (28) 
MeTui mgtui 

from which the result follows. H 



Remark. If there exists a vector ips invariant under all elements of R, it is 
invariant under all elements of T and satisfies \\ips\\ = IHV'slll = 1- 



5. Deterministic Results 

In this section, we derive some algebraic formulae and some uniform bounds 
for later purposes. Since there is no probabilistic issue involved here, we shall 
therefore simply denote Mj = M(u)j). We are concerned with the product 

>J?„ := Mi ■ ■ ■ M n . (29) 



5.1. Decomposition of the Mj 

With P\,Mj the spectral projection of Mj for the eigenvalue 1 we define 

<••.;: /'!.,<>•• Pr-=mm- (30) 

Note that {4>j \ips) = 1 so that Pj is a projection and, moreover, M*ipj = tpj . 
We introduce the following decomposition of Mj 

M, ■■ I', ■ (),M,(),. with Qj = l-Pj. (31) 

We denote the part of Mj in QjC d , by Mq j := QjMjQj. It easily follows 
from these definitions that 

PjP k - Pk, QjQk = Qj, (32) 
QjPk = o, PkQj = Pk - Pj = Qj - Qk- (33) 

Remark. If 1 is a simple eigenvalue, P\,Mj = Pj and (31) is a (partial) 
spectral decomposition of Mj. 
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Proposition 5.1. For any n, 

*n = \ips)(0 n \+M Ql ---M Qn , (34) 

where 

9 n = V„ + M£ n V„_i + • • • + M* Qn ■ ■ ■ M* 2 Vi (35) 
= M* ■ ■ ■ M 2 >! (36) 

and where (tl>s,0 n ) = 1. 

Proof. Inserting the decomposition (31) into (29), and using (32), (33), 
we have 

n 

Since Pj = \ips)(ipj\, this proves (34) and (35). From (33), we obtain for 
any j, k, 

M Q] M Qk = M Qj M k = QjMjMk. (37) 

Hence, tf „ = Pi Mi • • • M n + Q 1 M 1 ---M n = \^ s ) (M* ■ ■ • M|V>i | + 
M Ql ■ ■ ■ M Qn , which proves (36). □ 



5.2. Uniform Bounds 

The operators Mj, and hence the product are contractions on C d for 
the norm 1 1 1 • 1 1 1 . In order to study their asymptotic behaviour, we need some 
uniform bounds on the Pj, Qj, ■ ■ ■ Recall that \\ips\\ = 1- 

Proposition 5.2. Let Co be as in Corollary J^.l. Then, the following 
bounds hold 

(1) For any n E W, ||*„|| < C . 

(2) For any j E W, \\Pj\\ = \\^j\\ < C and \\Qj\\ < 1 + C„. 

(3) bup{||Mq, b M q , b _ i • ■■M Qn ||, n E N*, j fc £ N*} < C (l + Co). 
For any n G N*, ||0„|| < C^. 

Proof. It is based on Von Neumann's ergodic Theorem, which states that 

i 

Pi.Mj = lim - V Mj\ 

fe=0 
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The first two estimate easily follow, whereas the third makes use of (37) to 
get M Qjn M Qjn i ■ ■ ■ M Qji = Q jn M in M jn _ x ■ ■ ■ M h , so that 

W M Q in M Q in -i ■ ■ ■ M Qn II < WQinWCo < C (l + Co). 
Finally, (36) and the above estimates yield \\0 n \\ < CollV'ill < Co- ^ 



5.3. Asymptotic Behaviour 

We now turn to the study of the asymptotic behaviour of ^> n , starting with 
the simpler case of Ideal Repeated Interaction Quantum Systems. 
That means we assume 

M k =M, Vfc > 1. (38) 

If 1 is a simple eigenvalue of M, then P\ : m = \i>s)(i>\, for some tp s.t. 
(V#s) = 1, and 

*„ = M" = |Vs)(V'l+M2 (39) 

Further, if all other eigenvalues of M belong to the open unit disk, Mq 
converges exponentially fast to zero as n — > oo. 

Consequently, denoting by spr(iV) the spectral radius of N G Md(C), 

Lemma 5.1. If the RDO M belongs to M(e), 

*„ = |V5)(V'l+0(e" 7n ), (40) 
for all < 7 < 1 - spr(Mg). 

Two things are used above, the decay of Mq and the fact that 6 n — ip is 
constant, see (34). The following result shows that in general, if one knows 
a priori that the products of Mq^s in (34) goes to zero, converges if 
and only if P n = iV'sXV'nl) does. 

Proposition 5.3. Suppose that lim^oo sup{||MQ 3n ■ • • Mq h ||, jh G 

N*} = 0. Then 6 n converges if and only if ip n does. If they exist, these 
two limits coincide, and thus 

lim <J?„ = IV'sXV'ooI, 

n— >oo 

where V'oo = lim^oo ip n . Moreover, lips)^^] is a projection. 

In general, we cannot expect pointwise convergence of the 9 n , but we 
can consider an ergodic average of 9 n instead. This is natural in terms of 
dynamical systems, a fluctuating system does not converge. 
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The previous convergence results relies on the decay of the product 
of operators Mq. . Conditions ensuring this are rather strong. However, 
Theorem 6.1 below shows that in the random setting, a similar exponential 
decay holds under rather weaker assumptions. 

6. Random Framework 

6.1. Product of Random Matrices 

We now turn to the random setup in the framework of Section 3.2. For 
M(ui) an RRDO, with probability space (fi,.? 7 , p), we consider the RRDP 
on f2 N given by 

*„(S7) :=M{cj 1 )---M{cj n ), w G Q N *. 

We show that \l/ n has a decomposition into an exponentially decaying 
part and a fluctuating part. Let Pi(u>) denote the spectral projection of 
M(u) corresponding to the eigenvalue one (dim Pi (u ) > 1), and let Pi(co) 
be its adjoint operator. Define 

:= PiMVs, (41) 

and set 

P(u) = IVsXV'HI, Q(w) - 1 - P(u). (42) 

The vector tp(uj) is normalized as (ips, V'( w )) = 1- We decompose M(uj) as 

M(w) - P(w) + Q(w)M (w)Q(w) =: P(w) + M Q (w). (43) 

Taking into account this decomposition, we obtain (c.f. Proposition 5.1) 

*„(w) := • • • M K) = |Vs><fl n (w)| + • • • M Q (w„), (44) 

where 0„(o7) is the Markov process 

9 n (Zu) = M*(uo n ) ■ ■ ■ M*(lo2)iP(loi) (45) 

= VK) + M^(w„)V(w„-i) + • • • + Mq(w„) • • • M^(w2)^(wi), 

M*(u)j) being the adjoint operator of M{wj). We analyze the two parts in 
the r.h.s. of (44) separately. 

Theorem 6.1 (Decaying process). Let M(co) be a random reduced dy- 
namics operator. Suppose thatp(M(uj) G M(e)) > 0. Then there exist a set 
fil C £1 N and constants C\ a > s.t. P(fii) = 1 and s.t. /or any U G Oi, 
i/iere exists a random variable no{u>) s.t. for any n > n^(uj), 

||M Q ( Wl )-..M Q (a;„)|| <Ce~ an , (46) 

and E[e ana ] < oo. Moreover, E[M] G -M( B) . 
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Remarks. 1. The sole condition of M having an arbitrarily small, non- 
vanishing probability to be in M(e) suffices to guarantee the exponential 
decay of the product in (46) and that E[M] belongs to Mie\- 

2. Actually, E[M] G M( E ) is a consequence of spr(E[Mg]) < 1, which 
comes as a by product of the proof of Theorem 6.1. From the identities 

E[M] = |V's}W>]|+E[Mq], <E[V]|Vs> = 1, E[M q ]Vs = 0, (47) 

which do not correspond to a (partial) spectral decomposition of E[M], 
and this estimate, we get 

E[M]™ = \i/> s ){ E[ip] + E[Mq]*E[iP] + • ••E[Mq]*"- 1 E[^] \+E[M Q ]* n 

^oo - EIMq]*)"^^]! EE P 1MM] . (48) 

3. Our choice (41) makes 4>(u)) an eigenvector of M*{to). Other choices of 
(measurable) which are bounded in ui lead to different decompositions 
of M(lu), and can be useful as well. In particular, if M(u>) is a bistochastic 
matrix, can be chosen as an M*(w)-invariant vector which is indepen- 
dent of U). 



6.2. A Law of Large Numbers 

We now turn to the asymptotics of the Markov process (46). 

Theorem 6.2 (Fluctuating process). Let M(oj) be a random reduced 
dynamics operator s.t. that p(M(o>) G M(e)) > 0. There exists a set 0,2 C 
n ir s.t. P(f2 2 ) = 1 and, for all ZJ G VL 2 , 



1 N 

N lim co ^E^(-) = ^ ( 49 ) 



n=l 

where 

9 = Jim E[0 n ] = Pr, E[M] EM - P* MM] ^s- (50) 

Remarks: 1. The ergodic average limit of 8 n (ZJ) does not depend on the 
particular choice of ip(cu). This follows from the last equality in (50). 
2. The second equality in (50) stems from 

n-l 

E[» n ] = ^(E[Af Q ]fEM, (51) 

fe=0 

by independence, and which converges to P* E r M jE[^>] by (48). The third 
equality follows from (47). 
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3. Comments on the proof of these Theorems are provided below. 

Combining Theorems 6.1 and 6.2 we immediately get the following result. 

Theorem 6.3 (Ergodic theorem for RRDP). Let M(ui) be a random 
reduced dynamics operator. Suppose p(M(ui) 6 M.(e)) > 0. Then there 
exists a set ^3 C fi N s.t. P(^3) = 1 and, for all uj G ^3, 

1 N 

J im TvrZ! M ( t '>i)"- M K) = \4>s)(0\ =Pl,E[M]- (52) 

N—>oo 1\ ' 
n=l 

Remarks 1. If one can choose ip(u>) = ip to be independent of u>, then we 
have by (36) that 9 n (uJ) = if>, for all n,ZJ. Thus, from (44)-(46), we get 
the stronger result lim^oo M{loi) ■ ■ ■ M(u) n ) = \ips)(*P\, a.s., exponentially 
fast. 

2. This result can be viewed as a strong law of large numbers for the 
matrix valued process ^(w) = M{u\) ■ ■ ■ M(w n ). 

Comments: The existence of (ergodic) limits of products of random op- 
erators is known for a long time and under very general conditions, see 
e.g. 5 , 13 . However, the explicit value of the limit depends on the detailed 
properties of the set of random matrices considered. The point of our anal- 
ysis is thus the explicit determination of the limit (52) which is crucial for 
the applications to the dynamics of random repeated interaction quantum 
systems. 

The more difficult part of this task is to prove Theorem 6.1. The idea 
consists in identifying matrices in the product ^ n (uj) which are equal (or 
close) to a fixed matrix M that belongs to M(e)- Consecutive products 
of M give an exponential decay, whereas products of other matrices are 
uniformly bounded. Then one shows that the density of long strings of con- 
secutive M's in a typical sample is finite. Once this is done, a self-contained 
proof of Theorem 6.3, is not very hard to get 8 . 

On the other hand, given Theorem 6.1 and the existence result of 5 , 
we can deduce Theorem 6.3 as follows. Let us state the result of Beck and 
Schwarz in our setup. Let T denote the usual shift operator on f2 N defined 
by (Tw)j =u7j+i, j = 1,2,- ■•. 

Theorem 6.4 (Beck and Schwartz 5 ). Let M(u) be a random reduced 
dynamics operator on £1. Then there exists a matrix valued random variable 
L(ZJ) on £1 N , s.t. E[||L||] < 00, which satisfies almost surely 

L(U) = M(u}i)L(TU}), (53) 
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where T is the shift operator, and 

N 



Um _L £ M(wi) • • • MK) = L(ZJ). (54) 

N—>oc iv '-^ 

Further assuming the hypotheses of Theorem 6.1, and making use of the 
decomposition (44), we get that L can be written as 

L(uJ) = |tf s )(0(w)|, (55) 

for some random vector 8(uJ). Now, due to (53) and the fact that ips is 
invariant, dip) satisfies 

0(u7) = 9(TU) a.s. (56) 

The shift being ergodic, we deduce that 9 is constant a.s., so that 

0(u7) = E[0] a.s. (57) 

which, in turn, thanks to Proposition 5.2 and Lebesgue dominated conver- 
gence Theorem, allows to get from (46) 

E[0] = lim E[0„] =P* E[M] * 5 . (58) 
6.3. Limit in Law and Lyapunov Exponents 

We present here results for products "in reverse order" of the form <&„(w) := 
M(u) n ) ■ ■ ■ M(u)\), which have the same law as ty n (u>). They also yield infor- 
mation about the Lyapunov exponent of the process. The following results 
are standard, see e.g. 1 . The limits 

A*(u7) = lim ($„(w)*$„(o7)) 1 / 2 " and A*(cJ) = lim (* n ^)**„(cJ)) 1 / 2 " 

exist almost surely, the top Lyapunov exponent 71 (u;) of A$(ci7) coincides 
with that of A*(u7), it is constant a.s., and so is its multiplicity. It is in 
general difficult to prove that the multiplicity of 71 (w) is 1. 

Theorem 6.5. Suppose p(M(w) G M.{e)) > 0. Then there exist a > 0, a 
random vector 

?7oo(w)= lim ^(wi) + MQ(wi)^(w2) + --- + MQ(wi)---MQ(w n _i)^(w n ) 

and Q4 C fi N imi/i P(r^4) = 1 such that for any u G and n G N* 

$„(w)-|^s>(»?oo(w)|| <C^e- a ", /or some C-. (60) 
j4s a consequence, for any uJ G f^, 7i(w) is 0/ multiplicity one. 
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Comments: While the Theorems above on the convergence of asymptotic 
states give us the comfortable feeling provided by almost sure results, it is 
an important aspect of the theory to understand the fluctuations around 
the asymptotic state the system reaches almost surely. In our iid setup, the 
law of the product & n (u) of RRDO's coincides with the law of $„(uJ) which 
converges exponentially fast to \ips)(rj 00 (uj)\. Therefore, the fluctuations are 
encoded in the law of the random vector 7700 (tU). It turns out it is quite 
difficult, in general, to get informations about this law. There are partial 
results only about certain aspects of the law of such random vectors in 
case they are obtained by means of matrices belonging to some subgroups 
of Gld(R.) satisfying certain irreducibility conditions, see e.g. 10 However, 
these results do not apply to our RRDO's. 

6.4. Generalization 

A generalization of the analysis performed for observables acting on S only 
described above allows to establish the following corresponding results when 
instantaneous observables are considered. 

The asymptotics of the dynamics (20), in the random case, is encoded 
in the product 

M(wi) • • • M(w TO _ ; _i)./V(u; m _;, . . . , ui m+r ), 
where N : Q r+l+1 — ► M d (C) is given in assumption (R2). 

Theorem 6.6 (Ergodic limit of infinite operator product). 

Assume M(u>) is a RRDO and (R2) is satisfied. Suppose that p(M(w) G 
M(e)) 0- Then E[M] e M( E y Moreover, there exists a set fl 5 C fl N * of 
probability one s.t. for any U — (w«)„gn £ ^5, 

1 v 

lim -V]M(wi)-"M(w n )iV(w n+ i,...,w n+J+I . + i) = |V5><fl| E[JV], 

n—1 

where 9 = P^ E[M] ips- 

As in the previous Section, a density argument based on the cyclicity and 
separability of the reference vector ipQ allows to obtain from Theorem 6.6 
the asymptotic state for all normal initial states g on 971 given as Theorem 
3.3 
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